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Abstract 

The addition of non-renormalizable terms involving the Higgs fields to the MSSM (BMSSM) 
ameliorates the little hierarchy problem of the MSSM. We analyze in detail the two main cos- 
mological issues affected by the BMSSM: dark matter and baryogenesis. The regions for which 
the relic abundance of the LSP is consistent with WMAP and collider constraints are identified, 
showing that the bulk region and other previously excluded regions are now permitted. Requir- 
ing vacuum stability limits the allowed regions. Based on a two-loop finite temperature effective 
potential analysis, we show that the electroweak phase transition can be sufficiently first order in 
regions that for the MSSM are incompatible with the LEP Higgs mass bound, including parameter 
values of tan (3 <5, m^^ > mt, mg ^ TeV. 



I. INTRODUCTION 



The smallness of the quartic Higgs couphng in the framework of the minimal supersym- 
metric standard model (MSSM) poses a problem. The tree level bound on the Higgs mass 
is violated, and large enough loop corrections to satisfy the lower bound on the Higgs mass 
suggest that the stop sector has rather peculiar features: at least one of the stop mass 
eigenstates should be rather heavy and/or left-right-stop mixing should be substantial. 

The situation is different if the quartic Higgs couplings are affected by new physics. If 
the new physics appears at an energy scale that is somewhat higher than the electroweak 
breaking scale, then its effects can be parameterized by non-renormalizable terms. The 
leading non-renormalizable terms that modify the quartic couphngs are [H El [3l HJ |5], EJ [7] : 

Wbmssm = -^{HuHdf + -^Z{HuHdY J (1) 
where ^ is a SUSY-breaking spurion: 

Z = 6'^msusy (2) 

The first term in Eq. Q is supersymmetric, while the second breaks supersymmetry (SUSY). 
In the scalar potential, the following quartic terms are generated: 

2e,HMHiH^ + H\H,) + e^iH^H.f, (3) 

where 

The interplay between the Higgs sector, the stop sector, and the non-renormalizable 
(NR) operators has interesting consequences for the MSSM baryogenesis |8]. The window 
for MSSM baryogenesis is extended and, most important, can be made significantly more 
natural. In addition, these operators have implications for yet another cosmological issue, 
and that is dark matter [5]. In this work we present an extended analysis of the results for 
both the electroweak phase transition and the dark matter relic abundance in the BMSSM. 

One of the attractive features of the MSSM is the fact that the lightest R-parity-odd 
particle (LSP) is a natural candidate for being the dark matter particle. Progress in experi- 
mentally constraining the MSSM parameter space restricts, however, the regions where the 
dark matter is quantitatively accounted for to rather special regions of the MSSM: the focus 



point region, with surprisingly heavy sfermions; the funnel region, where the mass of the 
CP-odd neutral Higgs scalar is very close to twice the mass of the LSP; the co-annihilation 
region, where the mass of the scalar partner of the right-handed tau is very close to the mass 
of the LSP; and the bulk region, where the bino-LSP and the sleptons are light. 

The effects of the NR operators are potentially important for two of these four re- 
gions. First, these operators give rise to a new Higgs-Higgs-higgsino-higgsino interaction 
Lagrangian, 



which contributes to the annihilation process of two higgsinos to two Higgs particles. This 
effect is relevant when the dark matter particle has a significant component of higgsinos, as 
is the case in the focus point region. Second, as mentioned above, these operators modify 
the relation between the light Higgs mass and the stop masses. This effect can be important 
in the bulk region within models where the slepton and stop masses are related, such as the 
mSUGRA models. In this work, we will study these effects and assess their quantitative 
significance. 

The plan of this paper is as follows. In Section [IT] we present the BMSSM spectra of 
the Higgs, neutralino and chargino sectors, and the implications for the stop spectrum. In 
Section IIIII we describe the BMSSM modifications to the annihilation cross sections that 



are relevant to the dark matter relic abundance. In Section |IV| we analyze the implications 
of the BMSSM operators for dark matter, while in Section |V] we explore the parameter 
space where the electroweak phase transition (EWPT) is strongly first order, as required for 
successful baryogenesis. We summarize our conclusions in Section [VT} 
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II. THE SPECTRUM 



A. The Higgs sector and light stops 



We define the scalar Higgs components by 



u 



V2 
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h: 



4>2+Hur+iHu 
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(6) 



The VEVs of these Higgs fields are parameterized by 



tan/3 = 102/01 1, V = ^(0? + 0i)/2 ~ 174 GeV. 



(7) 



To leading order, the two charged and four neutral Higgs mass eigenstates are related to the 

interaction eigenstates via 
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where C/j = cos P, Sjs = sin /3, and similarly for a. Within the MSSM (without the 

operators), the angle a is given by (at tree level) 



S2a = 
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H 



mt 



S2I3- 



(9) 



If the ei,2 couplings are complex, then the four neutral mass eigenstates are related by a 4 x 4 
transformation matrix to the real and imaginary components of and H^. The analysis 
of dark matter is, however, unaffected at leading order, so we neglect such effects here. In 
the unitary gauge, the Goldstone fields and are set to zero. 



Taking mz-,mA and tan/5 as input parameters, we obtain the following ei-corrections to 
the Higgs spectrum: 
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The angle a is shifted from its MSSM value: 
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The possibility of a light Higgs scalar hidden in the LEP data is not excluded (see e.g. 
|lUj ) and, in fact, may have interesting implications for DM [llj. This scenario relies upon 
sizable mixing in the Higgs mass matrix, such that production via e'^e^ Zh is suppressed. 
Hence the mass splitting in the scalar Higgs sector must be rather small. A light Higgs 
sector is an interesting possibility also within the BMSSM, following the impact of the non- 



renormalizable operators on the mixing angles, captured to leading order in Eq. (11). An 



extreme example for this possibility was presented in [6] , in which the electroweak symmetry 
breaking vacuum is controlled by the non-renormalizable operators. 

In the present work, we limit our attention to the more conservative situation wherein the 
non-renormalizable operators can still be treated as perturbations in the usual electroweak 
breaking vacuum of the MSSM. We allow for significant splitting between the heavy and 
light Higgs mass eigenstates, in which case the constraint that arises from the LEP bound 
on the light Higgs mass can be written as follows: 



2(tree) 



ml ^ mr"""' + Stml + 5eml > (ll4 GeV)^ 



(12) 



where 



Stml 



+ 




47r2t>2 m? — m? 

tl 12 



Xt = At + n* cot p. 



(13) 



The Semi contribution relaxes the constraints on Stml in a significant way. In fact, 
St^l < is not excluded (in this respect we disagree with the conclusions of [7J). Thus, for 
ei < —0.05 and tan/? < 10, the two stop mass eigenstates can be as light as the top quark, 
or one could be as light as the direct experimental lower bound with the other only slightly 
heavier than the top. 

B. Neutralinos 



The neutralino mass matrix is given by 

^ Ml —mzSwC/3 mzSwS/3 ^ 

M2 mzcwcp -mzCwS/3 

-mzswcp mzcwcp — /i 

\ mzswsp -mzcwsp -fx 



+ 



w 



^0 ^ 







S2/3 



\o S2/3 4 y 

The transition to the neutralino mass basis is obtained with a unitary matrix Z: 



(14) 



= diag(m^^,mjY^,m^^,m^J Z. 



(15) 



The gaugino fraction in the LSP is defined as 



(16) 



while the higgsino fraction is given by 1 — i?^ = | Z13 p + | Z 
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C. Charginos 

The chargino mass matrix is given by 



M2 \/2mwSp 



1^0 ly 



(17) 



The transition to the neutralino mass basis is obtained with unitary matrices V and V: 



= [/^ diag(mg^,mgj V. 



(18) 
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III. ANNIHILATION CROSS SECTIONS 



The effects of the BMSSM operators on (co) annihilation cross sections are relevant to the 
dark matter issue when the LSP has a significant component of higgsino, namely 1 — ^ 1. 
These effects come in two ways. First, there is an indirect effect, due to the modification of 
the neutralino and chargino spectra. Since co-annihilation rates are very sensitive to mass 
splittings, this is often the more important effect. Second, there is a direct effect of the 
new ei-dependent couplings which modify the (co) annihilation processes that involve Higgs 
scalars as mediators and/or as final states. In this section, we focus on the latter effect. 



In Section IV we analyze numerically the DM relic abundance in the BMSSM, taking into 



account all effects. 

A. Single scalar 

We consider terms of the form 

C^(j)NN; C^(f)NC+, (f) = h, H, A; H' . (19) 

We denote the MSSM couplings by C° and define the modification that is induced by the 
ei terms as follows: 

C^ = C'^{l-6,^). (20) 

As concerns the emission of a single neutral scalar in the annihilation of two neutralino 
LSPs, we obtain the following 's: 

gM {Zi2-taiiewZn){saZi3 + CaZi4) 

^ _ 2y/2Xlv ^ cpCgZl^ + sps^Zl^ + 2g(a+/3) ^13^14 
gM {Zi2 - tan ewZu){-CaZri + SaZu) ' 
= v^^i^ X S2p{Zf, + Zf,)+4Z^,Z,, 

gM (Zi2 -tan 6*^/^11 )(s^Zi3 - c^Zu)' 

As concerns the emission of a single charged scalar in the co-annihilation of the neutralino 
LSP with the lightest chargino, we obtain the following ^^/^-'s: 

g ^ 2Ar^ ^ UujcpZu + s^Z.s) .241 

^^"^ gMtan(3 f/n^is - ^f/i2(^i2 + tan^^Zn) ' ^ ^ 

gM VuZ^. + j^VuiZu + tanewZiiY ^ ^ 
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The R and L sub-indices correspond to Pr and Pl which project a Dirac higgsino onto the 
lower and upper Weyl fermion, respectively. Thus, Ch-r is the coupling of the H~ N PrC^ 
term, while C^^-^, is the coupling of the H~ NP^C^ term. 

To get some intuition about the expected size of the correction to the MSSM annihila- 
tion cross section, we can estimate from the above expressions that the relative size of the 
correction, k: 

Thus, if there is a new physics threshold at around 5 TeV, and the higgsino component in 
the LSP is of order ten percent, then the correction to the annihilation cross section is of 
order ten percent. To make contact with the scalar spectrum, it is also useful to represent 



--o.el^lmfi^l . (27) 



Eq. (26) in terms of ei and /i, 

1/2 

/i ) vo.iy V 0.1 

As the new physics threshold scales like ~ /i/ei, it is difficult to envisage a phenomenolog- 
ically acceptable scenario exhibiting /i ~ m-w simultaneously with ei ~ 0.1. Hence some 



suppression is to be expected from the combination of mass and e factors in (27). This 
implies that the BMSSM modification to the relevant MSSM annihilation processes is at 
most of (9(10%) in generic cases. 



B. Two scalar s 



We consider terms of the form 



(28) 



where yah has mass dimension of —1. Keeping only potential s- and p-wave contributions, 
we obtain for the annihilation cross section into <^a'^h- 



ivihr+'jUy'atr+liyLr 



(29) 



where is the LSP mass. Cat = 2(1 + 6ab), S is the center-of-mass energy squared. 



(30) 
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and 

In terms of the Lagrangian parameters, we have 

fa, = -^^e (A*i;,) , = -^Irn (A*F,,) , 0,0, = hh, HH, AA, hH, R- 
= ^2:m (A*y,,) , yl^ = -l^ne{\\Yg,), c^M = hA,HA. (32) 
The dimensionless Yah couphngs are given by 

Yhh = sl^Zl^ + (?aZri — '^S2aZl4,Zi3, 
YhH = + S^Z^3 + 2S2a2'i4Zi3, 

= ~S^Z^4 — (?pZl^ — 2s2/32'i4Zi3, 
^/iH = ■52o( — ^14 + ^'^3) + 4C2o2'i4Zi3, 
^/lA = 2Sq,S/3Zi4 — 2CcjC/3Z^3 + 4S(q,_/3)Zi4Zi3, 

= — 2CaS/3Z^4 — 2SaC/3Z^3 — 4C(a-/3) ^'14^13 , 

= —'2s2pZi^Zi4. (33) 
Co-annihilation proceeds via terms of the form 

N'Mt + y:ti')ct, (34) 

with y'^, y° complex. The cross-section is given by an expression similar to Eq. (29 ), provided 



that the neutralino-chargino mass difference is neglected, taking Cab = 1, and making the 
substitution (y^)^ — ^ ll/'^P, (z/*)^ ^ We now have 

Vh-a = (KYh-ar T \iYh-al) , (35) 

yR-hR = V2Ci3Ui2{SaZi4 - C^Zis), 

YH-hL = V2Si3Vi2{SaZi4 - CaZis), 
Yh-HR = -V2CpUi2{CaZii + Sa^'is), 
Yh-HL = -V2Si3Vi2{CaZi4 + SaZis), 

Yh-ar = V2cpUi2{spZii + C/3Z13), 

Yh-al = V2spVi2{sf:iZi4 + cpZi^). (36) 



To get some intuition about the expected size of the correction to the MSSM annihilation 
cross section, let us consider the case of a higgsino LSP, R^^ ~ 0, wherein the couplings 



(33) and (36) are unsuppressed. In the MSSM, annihilation into (mostly transverse) gauge 

4 I I 

bosons gives (erf) ~ [12]. Regarding the BMSSM contribution (29) to annihilation 



into light Higgs boson pairs, we obtain (av) ~ 247rM^ * ^-^^^ estimate holds when there is no 
CP violation and the leading contribution is p-wave. Writing the modified cross section as 
(av) = (crt>)o(l + 5e), we find 6^ ~ 0.25 (f^)^- Note, however, that in the relevant scenario 
co-annihilations are important and so further numerical study is required to assess the full 
impact of the BMSSM. Below, we proceed to perform this study. 

IV. THE DARK MATTER RELIC DENSITY 

As deduced from the WMAP satellite measurement of the temperature anisotropics in 
the Cosmic Microwave Background, cold dark matter makes up approximately 23% of the 
energy of the Universe [I3]. The DM cosmological density is precisely measured to be 

(]^^^/i2 = 0.101 ±0.062 (37) 

at 68% CL. The accuracy is expected to be improved to the percent level by future mea- 
surements at Planck satellite |14] . 

We calculate the dark matter relic density in the presence of the ei couplings using a 
modified version of the code MicrOMEGAs [T3], where we implemented the BMSSM Higgs- 
Higgs-higgsino-higgsino couplings of Eq. (|5|. The leading ei^2-induced corrections to the 



MSSM Higgs spectrum, Eq. (10), were implemented using the code SuSpect [16j. 

The BMSSM framework, if relevant to the little hierarchy problem that arises from the 
lower bound on the Higgs mass, assumes a new physics scale at a few TeV. Since the new 
degrees of freedom at this scale are not specified, the effect of the new threshold on the 
running of parameters from a much higher scale cannot be rigorously taken into account. 
It therefore only makes sense to study the BMSSM effects in a framework specified at low 
energy. In order to demonstrate some of the most interesting consequences of the BMSSM 
operators for dark matter, we will employ two such sets of parameters: a model where all 
sfermion masses are correlated, and a model where the only light sfermions are the stops. 
The first model demonstrates how the so-called bulk region is re-opened, even for correlated 
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stop and slepton masses. The second model incorporates the interesting process of stop co- 
annihilation. For both models we focus our attention mainly on regions where the stops are 
light, since the main motivation for the BMSSM operators is to avoid a heavy stop (which is 
the cause of the little hierarchy problem) and, furthermore, this is the region that is relevant 
to BMSSM baryogenesis. Previous analysis in the context of the MSSM with a light stop 
were done in [HI [HI [H] . 

A. Correlated stop-slepton masses 

The most natural dark matter scenario within the MSSM framework could have been 
that of a light bino, annihilating to the standard model leptons via light slepton exchange. 
This scenario is known as the "bulk region" of the MSSM. However, in some of the most 
intensively studied MSSM scenarios, such as the mSUGRA [20] or cMSSM frameworks, the 
part of the bulk region that is allowed became smaller and smaller as the experimental lower 
bound on the Higgs mass became stronger. The generic reason for this is that a stronger 
lower bound on the Higgs mass requires a heavier stop which, in these frameworks, further 
implies heavy sleptons. One way to re-open the bulk region is to assume a framework where 
the stop and the slepton masses are not correlated. The BMSSM, however, re-opens the 
bulk region in a different way: the stop is not required to be heavy anymore. 

In order to understand these implications of the BMSSM framework and, in particular, in 
order to allow for a simple comparison with mSUGRA-like models, we investigate the follow- 
ing framework. The MSSM parameters that we use are those that would have corresponded 
to an mSUGRA model specified by the five parameters 

tan/3, mi/2, mo, Aq, sign(/i). (38) 

Thus, the correlations between the low energy MSSM parameters are the same as those 
that would hold in an mSUGRA framework. In other words, our low energy parameters are 
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expressed in terms of the parameters in (38 



) approximately as follows [2T] : 



m? ^ ml + 6m^/2, 
mj^ ^ ml + 0.5mi/2, 
"^L ^ ml + 0.15 ml/^, 

Ml ^ 0.4mi/2, 

M2 ^ 0.8mi/2, 

M3^3mi/2. (39) 

The values of fi^ and m\ depend on the soft breaking terms and on the electroweak breaking 
parameters in the standard way. Let us emphasize again that one should not think about 
this set of parameters as coming from an extended mSUGRA model, since the effects of the 
BMSSM physics at the few TeV scale on the running are not (and cannot) be taken into 
account. In addition, we have two extra BMSSM parameters: ei and €2- We focus essentially 
on the effects of ei. 

In practice, we make discrete choices of tan/5, Aq, sign(/i) and ei, and scan over uiq and 
fTT'i/2- We focus our attention on moderate values of mi/2 and rriQ because we are mainly 
interested in light sfermions and the bulk region. Fig. [l] displays the area, in the [mo, mi/2] 
plane, where the WMAP constraint is satisfied (between the solid red lines). We use Aq = 
GeV, /i > 0, and values of tan/?, ei and €2 as indicated in the various plots. The region 
below the dotted blue curve is excluded by the null searches for charginos and sleptons at 
LEP. The area to the left and above the orange curve is excluded because the stau is the 
LSP. The dotted black curves are contour lines for ruh. 

We would like to emphasize several points regarding the effects that are demonstrated 
by Fig. [1} 

1. The most significant effect of the BMSSM operators is their impact on the Higgs 
boson mass. Within the MSSM with mSUGRA-like correlations, the bound on the 
Higgs mass strongly constrains mi/2. In contrast, in the presence of ei = (9(— 0.05), 
the full region for which the correct value of the relic abundance is obtained is allowed. 

2. Another significant effect is the appearance (for negative enough ei values) of a new re- 
gion fulfilling the DM constraint. This is the '/i-pole' region in which rrih ~ 2m^o, and 
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FIG. 1: Regions in the [mo, '^'1/2] plane in which the WMAP constraint is fulfihed (between the 
sohd red hnes). We use ^0 = GeV, /U > 0, and values of tan /?, ei and £2 as indicated in the various 
panels. The region below the dotted blue curve is excluded by the null searches for charginos and 
sleptons at LEP. The area to the left and above the orange curve is excluded because the stau is 
the LSP. The dash-dotted black curves are contour lines for m/j with values in GeV as indicated. 
The dotted purple line in the middle right panel denotes the appearance of a remote vacuum in 
the scalar potential. Above and to the left of this line, the electroweak vacuum is metastable. 
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the s-channel Higgs exchange is nearly resonant, allowing the neutralinos to annihilate 
efficiently 



3. For the case discussed in this subsection, the co-annihilation region with the stops is 
not present as the mass difference with the LSP is always too large. 

4. In the mo region that we are considering here, the impact of the BMSSM operators 
on the mass of the neutralino LSP and on the region fulfilling the WMAP constraint 
is rather limited. The reason is that in the bulk region the LSP is mostly bino-like, 
while the BMSSM operators affect the higgsino parameters. 

The BMSSM operators may destabilize the scalar potential. If 4|ei| > €2, the effective 
quartic coupling along one of the D-fiat directions is negative, causing a remote vacuum 
to form in the presence of which the electroweak vacuum could become metastable. When 
considering values of ei > —0.1, vacuum stability is ensured provided that the following 
condition is fulfilled [23j: 

The stability criterion as written here applies when > 3m z- The fact that this criterion 
involves electroweak scale parameters, and only ratios of non-renormalizable operators, is 
reminiscent of the supersymmetric origins of ei. In regions of the parameter space where 



Eq. (40) is violated, the electroweak vacuum is metastable. For large tan/3, this occurs 
throughout an important fraction of the [mo, mi/2] plane, depicted in the middle right panel 
of Fig. [l]by the area above and to the left of the dotted purple line. 



It is important to stress that Eq. (40 ) represents an analytical tree level approximation and 
as such, in the current framework where quantum corrections are sizeable, it is conservative. 
In order to determine whether the lifetime of the electroweak vacuum is long enough, the rate 
of quantum tunneling into the remote vacuum should be compared to the Hubble rate. For 
illustration, consider the two parameter points. Pi : [mo = 200 GeV, mi/2 = 150 GeV] and 
P2 : [mo = 100 GeV, mi/2 = 250 GeV], with tan (3 = 10, ei = -0.1, €2 = 0, as in the middle 
right panel of Fig. [Tj Point Pi corresponds to a stable electroweak vacuum configuration. 



in agreement with the criterion (40). Point P2, on the other hand, violates Eq. (40) by 
~ 15%. However, P2 involves significant quantum corrections due to stops (m^-^ = 415 GeV, 
= 590 GeV, At = —480 GeV). Calculating the tunneling action, we find that P2 is in fact 
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long-lived enough to provide an acceptable model point. Similar calculations reveal that the 



rest of the middle right panel is also long-lived enough (at least marginally), and so Eq. (40) 
can not be used to exclude regions in the mSUGRA-like parameter space considered in this 



section. We will return to the issue of vacuum stability in Section IV B There, quantum 
corrections due to stops will be held moderate and fixed, resulting with precise application 



of Eq. (40) to exclude significant portions of the parameter space. 

As concerns precision electroweak data and low energy processes, it is important to realize 
that the new physics that generates the non-renormalizable operators can directly modify 
the constraints that come from these measurements. Ignoring this point, it is still possible 
to identify regions in the parameter space favored by the WMAP data which satisfy all 
such low energy constraints. The relevance of the BMSSM lies in the fact that constraints 
involving the Higgs are decoupled from constraints involving the stop sector. In particular, 
stops that are neither heavy nor mixed are acceptable, as demonstrated by our choice of 
Aq = throughout the current section. In contrast, in the case of the MSSM, satisfying the 
Higgs mass bound as well as electroweak precision data necessitates large values for Aq 



B. Light stops, heavy sleptons 

The aim of this section is to further expose various implications of the BMSSM for 
the DM relic abundance, putting special emphasis on parameter regions compatible with a 
strong first order electroweak phase transition, as required for baryogenesis. In particular, 
we are interested in the scenario of light, unmixed stops. As mentioned above, LSP co- 
annihilation with stops is not a viable possibility if the low energy soft supersymmetry 
breaking parameters obey relations similar to those that would follow from mSUGRA-like 
theory, unless stop mixing is very large. To explore this possibility in the BMSSM, we employ 
a set of low energy parameters that is different from the previous subsection. Explicitly, in 
addition to the BMSSM parameters, we consider the following set of parameters: 

M2, /i, tan/3, Xt, mu, mQ, mj, uia, (41) 

where is a common mass for the sleptons, the first and second generation squarks, and 
Ir. We further use Mi = | tan^ ^vk M2 ~ | M2. 

To demonstrate our main points, we fix the values of all but two parameters as follows: 
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= or -0.1, ea = or +0.05, tan/3 = 3 or 10, Xt = 0, mu = 210 GeV, tuq = 400 GeV, 
mj = rriA = 500 GeV. This scenario gives rise to relatively light stops: 

rrii^ < 150 GeV, 370 GeV < m^-^ < 400 GeV. (42) 

We scan over the remaining two parameters, M2 and /i. 

In the prescribed framework, one can identify four regions in which the WMAP constraint 
is fulfilled: 

- The 'Z-pole' region in which the LSP is very light, m^o ~ ^Mz ~ 45 GeV, and the 
s-channel Z exchange is nearly resonant. This region is not ruled out only in scenarios 
where the mass splitting between Mi and M2 at the electroweak scale is very large. 

- The '/i-pole' region in which the LSP is rather light, m^o ~ ^M^, and the s-channel h 
exchange is nearly resonant, allowing the neutralinos to annihilate efficiently [22]. 

- The 'mixed region' in which the LSP is a higgsino-bino mixture [12], M2 ~ 2fi, which 
enhances (but not too much) its annihilation cross-sections into final states containing 
gauge and/or Higgs bosons: XiXi ~^ W~^W~ , ZZ, Zh and hh. 

- The 'stop co-annihilation' region, in which the LSP is almost degenerate in mass with 
the lightest stop (ti). Such a scenario leads to an enhanced annihilation of sparticles 
since the Xi ~ ti co-annihilation cross-section [2Sl ES] is much larger than that of the 
LSP. 

Fig. [2] displays the areas, in the [M2, fi] plane, in which the WMAP constraint is satisfied 
(between the solid red lines). The region of large fi and M2, to the right of the orange dashed 
line, is excluded since the lightest stop becomes the LSP. The region of small fi and/or M2, 
below and to the left of the blue dotted line, is excluded by the null search for charginos at 
LEP2 [27]. 

Let us first consider the two upper panels, where Ci = €2 = 0. The various regions 
described above that are consistent with the DM constraints can be identified in this figure. 
The region at M2 ~ rriz ~ rrih corresponds to the s-channel exchange of an almost on-shell 
Higgs or Z boson. Note that when 2m^o is too close to the Higgs or Z mass pole, the 
LSP annihilation is too efficient and leads to a much too small ^dm h!^- In any case, for the 
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FIG. 2: Regions of the [M2, /i] plane in which the WMAP constraint is fulfihed (between the 
sohd red hnes). We use the parameter set described in the text, and values of tan/3, ei and £2 
as indicated in the various panels. The region above and to the right of the orange dashed curve 
corresponds to a stop LSP. The region to the left and below the blue dotted curve is excluded 
by direct chargino searches. Most of the region above the horizontal dashed line is excluded by 
vacuum instability. 
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Higgs mass values obtained here, nih ~ 85(98) GeV for tan/? = 3(10), this region is aheady 
excluded by the negative searches for chargino pairs at LEP2. 

The region close to M2 ~ /i ~ 200 GeV corresponds to the LSP being a bino-higgsino 
mixture with sizeable couplings to W, Z and Higgs bosons, allowing for reasonably large 
rates for neutralino annihilation into XiXi ~^ W~^W~ , ZZ, hZ and hh final states. Above 
and below the band, the LSP couplings to the various final states are either too strong or 
too weak to generate the relevant relic density. 

Finally, for larger fi values, in the region close to the orange curve, the mass of the 
lightest neutralino approaches the mass of the lightest stop leading to an enhanced co- 
annihilation cross-section: Xi^i ~^ b, gt {r^ 90%). Also, to a lesser extent (~ 5%), 
the annihilation cross-section of the stop NLSP contributes to the total cross-section by the 
process ti ii ^ g g. 

Next we consider the ei = —0.1 case (the two middle panels). The features of the 
DM allowed regions are similar to the previous case. The main difference comes from the 
important enhancement of the Higgs mass due to the presence of the BMSSM operators. In 
this case it is possible to disentangle the Z and the h peaks, since the Higgs-related peak 
moves to higher M2 values, due to the increase of the Higgs mass: rrih = 122(150) GeV for 
tan/3 = 10(3). Furthermore, the latter peak is no longer excluded by chargino searches. 

For large values of the /i parameter, the BMSSM operators destabilize the scalar potential. 
The regions above the horizontal dashed lines in Fig. [2] exhibit a metastable electroweak 



vacuum, following from the violation of Eq. (40) which requires fi < rriA- Similarly to the 
situation in the previous section, points right above the approximate analytical stability 
line may still be acceptable, as the stability constraint is somewhat alleviated by quantum 
corrections. However, with stop parameters as specified, computing the tunneling action 



for the middle and lower right panels reveals that the limit dictated by Eq. (40) is in fact 
accurate to better than 15%. Hence a significant portion of the parameter space in the 
middle and lower right panels is indeed excluded by stability considerations. 

The role of €2 can be seen from the lower two panels. It alleviates the stability constraint 
and slightly increases the Higgs mass. The e2-related effect on the Higgs mass is suppressed 
by tan^ (3 and consequently it is much more pronounced in the tan (3 = 3 case than in the 
tan/5 = 10 case. 



18 



V. THE EWPT IN THE BMSSM 



To study the electroweak phase transition, we analyze the finite temperature effective 
potential at two-loop order for the light scalar field. Detailed analyses for the case of 
the MSSM have been performed in Refs. [281 EHJ [301 El]- Here we focus on the region 
of parameter space of large niA, for which there is a single light Higgs field at the phase 
transition. We review the main effects found in Ref. [8J and study in more detail the interplay 
between the relevant parameters. We write the effective potential at finite temperature as 
follows: 

7 - 51n 



T 

The leading contributions to the 7-term and the E'-term arise at one-loop. The 7-term 
is further corrected at two- loops by the so-called Dsv and Dss terms [28]. The 5-term 
contains the dominant two- loop corrections in the MSSM, which arise from the Dssv and 
Dsss diagrams that contribute with a logarithmic dependence [2S]. For the numerical 
analysis we keep, in addition to the SM contributions, all contributions associated with the 
light right-handed stop. When considering a light left-handed stop, the relevant, dominant 
two-loop (order g'^h'f and hf) corrections are added as well. 

The strength of the phase transition is determined by the order parameter which, at the 
critical temperature of the phase transition T^, is given by 

0(T,) E 1{E^ 25 V^' E \B\ 

+ T7 + ^ 1 + 1^ • (44) 



T, 2A 2 V A2 A y X \ Ey 
We now discuss the various relevant parameters {E, B, A) and their dependence on the MSSM 
parameters. 

The cubic term arises in the regime for which the high temperature expansion of the 
bosonic one-loop contribution to the effective potential is valid. We write 

E = EsM + -Emssm- (45) 

The SM contribution is given by 

Esu = 2!!|±!4. (46) 

The dominant MSSM contribution comes from the stops (we provisionally neglect the 1^ — 1^ 
mixing) : 

^^ = -^4^{!4ST) + rnl{T)). (47) 
19 



The finite temperature stop masses are given by 

1 



(48) 



where we neglect contributions from the f/(l)y gauge bosons in the field dependent terms. 
To maximize the value of -Emssm, one would like to take negative values of mg and mfj, 



rrir 



m 



u 



that would cancel the thermal masses, yielding a purely cubic form in [29]: 



-^MSSM 



(49) 



(50) 



Eq. (50) illustrates the effect of the stops and gives, to leading order, what would be 
an upper bound on the strength of the phase transition for the MSSM. It is impossible, 



however, to make the selection (49) simultaneously for both stops, due to constraints from 
the p parameter and the experimental bound on the sbottom mass (through its dependence 
on mg). 

The two-loop stop contribution to the finite temperature potential, which can increase 
0(Tc) via its effect on B, is given by 



6V 



9l{Nl - 1)T^ 
167r2 



(T) log 



2m,- JT) 
3T 



+ m?jT) log 



2m,- JT) 
3T 



The effect is maximal when the limit of Eq. (49) is realized: 



B. 



stops 



(51) 



(52) 



w 



Other logarithmic terms tend to diminish the value of B [28j. Note, however, that the net 



two- loop contributions to Eq. (44) is the same when considering the maximal contributions 
from both stops or from a single one of them. 

The strength of the phase transition is further affected by the quartic coupling A. At zero 
temperature, it is related to the Higgs mass via ml = Aefff^, where 



— m^Cg^ — 16f ei cot fj + -^^mf.hf. m — 2 ■ 



m 



(53) 



Eq. (53) is valid for large itla and large tan/5. It includes the leading one- loop corrections. 
The €2 dependence is dropped. Adding in the leading finite temperature correction, we have 

3 



A(T) = Aeff-f ^/^^Jln2. 



(54) 
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FIG. 3: Contour of (pc/Tc = 1 (solid black curve), and contours of rrih (dashed curves) for tan (3 = 5, 
and ei = —0.065 (left panel) or —0.08 (right panel). Also shown is the constraint arising from the 
lower bound on m^^ (solid blue line). 



Thus, we can estimate. 



<h _ v'^jEsM + -^mssm) 



(55) 



We now consider three specific cases. We use different values of the non-renormaUzable 
contribution (ei) and the loop contribution (m^^m^^), but in such a way that rrih is fixed at 
the experimental lower bound. 

(a) m^f^"" = 114 GeV: 



To keep rrih = 114 GeV (and by that minimize the denominator of Eq. (55)), we need the 
loop correction of the Higgs boson mass to vanish, namely 



(56) 



This constrains the possible values of m'^^{T) and m'^^{T), and gives the largest possible 
contribution to -Emssm and to B, thus maximizing the strength of the phase transition. 
Smaller values of nii^nii^ are not allowed due to the experimental bound on nih. Larger 
values are allowed but lead to a decrease in the numerator and an increase in the denominator 



in Eq. (55). 



The left panel of Fig. [s] shows the contour of (pc/Tc = 1 and the contour of nih = 
114 in the mi^-mi^ plane. We also show the curve of m^^ > 175 GeV, corresponding to 



the experimental lower bound on m^^. We use tan/3 = 5 and ei 



-0.065, so that we 



obtain w}^'^'^ = 114 GeV. The rrih curve then corresponds to Eq. (56). The plot shows the 
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approximate symmetry in this case between m^^ and rrif:^, which is broken only by effects 
arising from D-terms. The allowed region is to the left of the (pc/Tc = 1 curve and above 
the rrLfi = 114 GeV curve. 

(b) m)^'''' > 114 GeV: 

Here we take values of ei such that m)^'^'^ > 114 GeV. In this case, m^^mi^ < mf is allowed. 
Such values further increase -Emssm with respect to the values presented in case (a), and the 
allowed region for a strong enough phase transition is increased, although here the sbottom 
mass constraint eliminates a large portion of this region. The right panel of Fig. [3] shows 
the (pc/Tc = 1 the rrih = 114 GeV contours for this case. 

(c) m)^'''' < 114 GeV: 

This is reminiscent of the usual MSSM results for the electroweak phase transition, in which 
one requires m^^ to be large enough to satisfy the m'j^'^ value. The effect on -Emssm is 
to effectively screen the contribution from the left-handed stop, thus reducing the cubic 
term to be smaller than half its maximum value and, furthermore, reducing the value of 
B. Variations in the value of mg produce only small variations in the value of rriu, as the 
main effect of increasing mq is to increase the one-loop value of the Higgs boson mass, thus 
reducing the strength of the phase transition. To compensate for that, a slightly smaller 
value of m^^ must be taken to increase -Emssm- 

In the region of the parameter space where the EWPT is strong enough, a minimum 
where colour is broken might develop [22] • If the temperature where this minimum becomes 
as low as the potential at the origin, Tj/ , is higher than the critical temperature for the 
EWPT, Tf, then the Universe is likely to end up in the colour breaking minimum. Thus, 
this region is excluded [32]. When the two stops are light enough, is safely higher than 
TJ/ . However, when we consider higher and higher values of ttl^^, and correspondingly, to 
keep the EWPT strong enough, lower and lower values of ^t^, the closer do and get 
to each other, until, at some critical values of (mS ,mS ), we reach T't = TY . 

Up to this point, we fixed the values of ei and tan/3 and obtained the allowed regions in 
the {m^^jnif:^] plane. We learned that the presence of ei = O{—0.1) opens new regions in 
this plane where both the rrih constraint and the EWPT constraint are satisfied. We now 
turn our attention to the dependence on the other relevant parameters. More concretely, 
we fix rrij^ = 200 GeV, which is close to the lowest value allowed by the 6p parameter, and 
obtain the dependence of (pc/Tc and of on tan/3,mr , and ei. 
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Our results are presented in Figs. |4]as contours of (j)c/Tc and of rrih in various parameter 
planes. In the top left panel we make discrete choices of ei, and present the results in the 
[mj-^,tan/?] plane. For each value of ei, the allowed region is to the right of the = 114 
GeV and to the left of the (pc/Tc = 1 curve. In the top right panel we make discrete choices 
of and present the results in the [ei,tan/5] plane. For each value of mj^, the allowed 
region is to the left of the rrih = 114 GeV and to the right of the (pc/Tc = 1 curve. In the 
lower panel we make discrete choices of tan/3 and present the results in the [ei, rrif:^] plane. 
For each value of tan/3, the allowed region is to the left of the rrih = 114 GeV and below 
and to the right of the (pc/Tc = 1 curve. We learn that small values of tan (3 are allowed and, 
furthermore, values of m^-^ > rrit can simultaneously satisfy the requirement for baryogenesis 
and the Higgs boson mass bound. 

To summarize the results depicted in Figs. [3] and |4| the BMSSM framework introduces 
the following significant differences with respect to the standard MSSM case: 

1. The Higgs boson mass constraint decouples from the value of mg in some regions of 
parameter space. Consequently, a very light left-handed stop can contribute signif- 
icantly to the cubic term in the effective potential. Explicitly, values of m^-^ ^175 
GeV {niQ ~ 0) are now allowed, compared to the values of tuq ~ 2 — 3 TeV, that are 
required in the MSSM. 

2. The lowest allowed value for m^^ arises from electroweak precision measurements and 
from direct searches. For the MSSM, these constraints are superseded by requiring 
nih > 114GeV. 

3. In the MSSM, a large value of left-right mixing in the stop sector is strongly favored. 
This is not necessary within the BMSSM. The actual consequences of left-right mixing 
in the stop sector are quite similar between the BMSSM and the MSSM: The Higgs 
boson mass is increased, which, in turn, weakens the phase transition (although this 
can be compensated by changing ei). Moreover, the contributions to the cubic term 
in the potential from the stops are screened, again weakening the phase transition. 

4. Larger values of m^^ are allowed. Actually, the experimental bounds on the masses of 
the left handed squarks determine the largest possible mass of the right-handed stop. 

5. Smaller values of tan (3 are now feasible. 
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FIG. 4: The dependencies of the Higgs mass and the EWPT on tan /3, ni^^, and ei. We fix 



mr 



200 GeV. Solid hnes indicate contours of (t)c/Tc = 1 and dashed hnes contours of rrih = 114 



GeV. Each same-colour pair of (solid and dashed) lines corresponds to the same set of parameters. 



We learn that the greater freedom in the values of the different parameters that affect 
the Higgs boson mass - tan/?, mi^, m^^ - implies that the allowed regions are similar to 
those identified in the early analyses of the EWPT |30j . 

Two final comments are in order. First, note that the inclusion of non-renormalizable 
operators of dimension six could modify the results. In particular, it is clear that for large 
rriA, such terms can induce additional contributions like those identified in Ref. [33]. Vacuum 
stability considerations must be carefully taken into account in this case. 

Second, as pointed out in Refs. |5l [8], the non-renormalizable operators can induce new 
sources of CP violation in the scalar Higgs sector, that would modify the production mech- 
anism of the baryon asymmetry of the Universe. These new sources could provide a relief to 
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the tension between the large phases required to produce the BAU and the contributions to 
the electric dipole moments. Additionally, even in the absence of CP violation in the scalar 
Higgs sector, modifications to the EDMs arise from the new interaction terms between the 
Higgs bosons and the higgsinos [4j. These issues will be further studied elsewhere. 

VI. CONCLUSIONS 

The main motivation to add non-renormalizable operators to the MSSM Higgs sector 
is to reduce the fine-tuning that is required by the lower bound on the Higgs mass. We 
have shown that, in addition, these operators have implications for cosmology that are 
very welcome: Regions of the super symmetric parameter space that are favored by the 
dark matter constraints and by the requirement for a strong first-order electroweak phase 
transition, but are excluded within the MSSM, become viable within the BMSSM. 

As concerns dark matter, a particularly important feature of the BMSSM is that the 
bulk region, for which the LSP is mostly bino-like and the sfermions are relatively light, 
can provide the adequate contribution to the energy density of the Universe while still 
satisfying the collider constraints on the Higgs boson as well as on the supersymmetric 
particle spectrum. Light stops co-annihilating with the LSP could have been active players 
in driving the dark matter relic density to its present value. It is also possible that nearly 
resonant LSP annihilation proceeded through exchange of the lightest Higgs particle itself. 

If light stops are indeed found in upcoming experiments, large BMSSM corrections will 
be implied. In this scenario, parameter regions where fi is large (exhibiting some heavy 
neutralinos and charginos) will be significantly constrained by the requirement of vacuum 
stability. 

As concerns the electroweak phase transition, the BMSSM has a dramatic effect when 
determining the range of parameters for which the phase transition is sufficiently strong to 
suppress sphaleron transitions in the broken phase. The fact that large stop-related radiative 
corrections to the Higgs mass are not required, allows light stop degrees of freedom to affect 
the dynamics of the phase transition by enhancing their contributions to the magnitude of 
the order parameter at the critical temperature. 
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Note added: 

Upon completion of this work, a related paper [31] has appeared. Where the two analyses 
overlap, we confirm their results. 
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